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We use first-principles calculations to extract two essential microscopic parameters, the charge- 
transfer energy and the inter-cell oxygen-oxygen hopping, which correlate with the maximum super- 
conducting transition temperature T c , max across the cuprates. We explore the superconducting state 
in the three-band model of the copper-oxygen planes using cluster Dynamical Mean-Field Theory 
with an exact diagonalization impurity solver. We find that variations in the charge-transfer energy 
largely accounts for the trend in T Cjmax across the cuprate families. 



PACS numbers: 0.0 

Despite an immense body of theoretical and exper- 
imental work, we have limited microscopic insights of 
which materials-specific parameters govern the trends 
in the maximum transition temperature Tcmax across 
the copper oxide superconductors. Structurally, all the 
cuprate families have in common C11O2 planes which sup- 
port superconductivity. They are described by the chemi- 
cal formula YS , n -i(Cu02) n , where n CuC>2 planes are in- 
terleaved with n—l spacer layers S to form a multi-layer. 
These multi-layers are then stacked along the c-axis, sep- 
arated by a different spacer layer X. Empirically, it is 
known that T C)lnax is strongly materials-dependent, rang- 
ing from 40 K in La2Cu04 to 138 K in HgBa2Ca2Cu30s. 
Additionally, Tc.max can be tuned both as a function of 
doping and the number n of Cu02 planes. 

Studies linking the known empirical trends to micro- 
scopies have generally established that the properties of 
the apical atoms (O, F or CI, depending on the cuprate 
family) are the relevant materials-dependent parameters. 
However, conclusions vary regarding their effects on elec- 
tronic properties, especially in multi-layer cuprates where 
not all Cu02 have apical atoms. Early theoretical work 
by Ohta, et. al., found correlations between T c and the 
Madelung potential of the apical oxygen, arguing that 
the apical potential controls the stability of the Zhang- 
Rice singlets [1]. They conclude that ^Cu-O' * ne distance 
between the Cu and apical O, is uncorrelated with su- 
perconductivity. In a more recent DFT study, Pavarini, 
et. at, argue that tunes between the single-layer 

cuprate families, affecting the electronic structure pri- 
marily via the one-electron part of the Hamiltonian [2]. 
Moving the apical oxygens away from the copper oxide 
plane allows stronger coupling of in-plane O 2p orbitals 
to the Cu 4s, enhancing the strength of longer ranged 
hoppings. This effect is characterized by the increase 
of a range parameter r ~ tf /t, describing the relative 
strength of the next-nearest neighbor hopping if to near- 
est neighbor hopping t in a one-band model. They find 
that materials with larger r have larger T Cjmax . Many- 
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FIG. 1. Parameters of the three-band p-d model for the CUO2 
planes in the cuprate superconductors. We show the two 
shortest-ranged oxygen-oxygen hoppings t vv and t pp i , and the 
on-site energies ed and e p . 

body corrections to t' were included by Yin, et. al. [3] . 

The development of cluster Dynamical Mean-Field 
Theory (c-DMFT) combined with first-principles calcu- 
lations [4] has advanced our qualitative and quantitative 
understanding of the cuprates [5, 6]. A satisfactory de- 
scription of these materials at intermediate temperature 
and energy scales has been achieved, and the consensus 
is that the cuprates lie in the regime of intermediate cor- 
relation strength [7-9] near the Zaanen-Sawatzky- Allen 
(ZSA) boundary [10]. 

In this Letter, we address the origin of the variation of 
the experimental T c . max across the cuprates using recent 
advances in electronic structure methods. We carry out 
first-principles calculations of the hole-doped cuprates 
and downfold to the 3-band p-d model to extract band- 
structure parameters. We correlate T Cjmax against the 
extracted parameters and organize the cuprate families 
by the magnitude of their charge-transfer energy. Us- 
ing c-DMFT, we explore the superconducting state as a 
function of the relevant microscopic parameters, result- 
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FIG. 2. In single-layer cuprates, increasing the apical oxygen 
distance reduces the charge-transfer energy. 



ing in suggestions for possible improvements in materials 
design to reach higher critical temperatures. 

Effective low-energy hamiltonians containing the min- 
imal set of bands are important tools for understanding 
chemical trends. We use the Wien2K code [11] to perform 
Linearized Augmented Plane Wave (LAPW) calculations 
on all major copper oxide families, and then extract 
model hamiltonian parameters by downfolding [12] to or- 
bitals constructed in the manner described in Ref. 13. 
In this work, we choose to downfold to a 3-band hamil- 
tonian describing the in-plane Cu-3d x "-y" and 0-2p or- 
bitals (Fig. 1). We believe four parameters capture the 
essential physics: the charge-transfer energy — e p be- 
tween the Cu and O atoms, the direct Cu-0 hopping 
t p d and the two shortest-ranged 0-0 hoppings t pp , and 
tpp/ . The extracted values are tabulated in Table I of the 
Supplementary Material. 

We find that only two parameters, — e p and t pp >, 
vary significantly across the cuprates. Although not cru- 
cial for our subsequent work, one would like to have a 
simple structural explanation for these trends. For the 
single-layer cuprates, the variation can be directly con- 
nected to dcu_°o (also tabulated in Table I). As we bring 
the negatively-charged apical oxygen towards the CuO 
plane, the resulting electrostatic repulsion suppresses the 
hopping t pp i , since t pp ' describes transitions of electrons 
past the Cu site, and provides justification for fact that 
tppi is smaller than t pp [14]. This mechanism for the de- 
pendence of hoppings on dg^Q has been pointed out in 
Ref. 2 for one-band models. However, we show in Fig. 2 
that the electrostatic repulsion simultaneously increases 
e<i ^ £p by rendering it costly to place an electron on the 
Cu site. These simple structural trends are less clear for 
multi-layer cuprates, where additional variables such as 
the inter-layer distance introduce additional complexity. 

Having identified the two relevant parameters, we plot 
7c,max against these quantities in Fig. 3a and Fig. 3b to 
identify possible correlations. Beginning with La2Cu04 
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FIG. 3. Correlations of T Cjmax in the copper oxides with the 
microscopic parameters of the three-band model Hamiltonian 
with (a) the charge-transfer energy e^ — e p (b) the next-nearest 
neighbor oxygen-oxygen hopping t pp i (c) the effective one- 
band range parameter r ~ t' /t. The trend of the dependence 
of the one-band range parameter agrees with Ref. 2. 



(LSCO), the limiting case among the cuprates since it has 
the largest e ( i — e p as well as the smallest t pp i , the figures 



show that both (i) decreasing — e p and (ii) increasing 
tp p i correlates with a enhanced T c max . To relate our 
results with the one-band model, we further integrate 
out the oxygen orbitals and compute the range parameter 
r ~ t'/t. Our plot of Tcmax against r in Fig. 3c agrees 
with the correlation found in Ref. [2]. 
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In order to clarify how the identified microscopic pa- 
rameters control T c ,max, we use c-DMFT in the cellu- 
lar form [4] with a 2 x 2 cluster of impurities to solve 
the downfolded three-band model. The non-local self- 
energy in c-DMFT captures the short-ranged correlations 
which are crucial to describe <i-wave superconductivity. 
Since the fermionic minus sign problem prevents finite- 
temperature solvers from accessing the low-temperature 
superconducting regime, we use an exact diagonalization 
(ED) impurity solver [15]. In this work, we extend previ- 
ous c-DMFT calculations of the one-band model [16, 17] 
to the three-band model, with realistic parameters ob- 
tained from first-principles calculations. The refinement 
captures the admixture of the Cu and O character near 
the Fermi level via a bath representing both the Cu and 
O degrees of freedom in the DMFT self-consistency con- 
dition. 

The three-band hamiltonian we treat with c-DMFT is 
as follows: 

H= tffc\ a<T c j/3lJ +^2e a n la<J + U dd J2 n id^n idi 

where i,j run over the in-plane CuC>2 unit cells, a, j3 
label the orbitals p x , p y and d x 2_ y 2, and a is the elec- 
tron spin. The hoppings tfr and onsite energies e a are 
those sketched in Fig. 1, except for the d-orbital onsite 
energy, where we subtract out a doping- and material- 
independent double-counting correction E^c to remove 
the correlations already captured within LDA [18]. The 
d x 2_ y 2 on-site Coulomb repulsion is U dd — 8 eV, although 
the results are not sensitive to the exact value as long as 
it is larger than the value of the charge-transfer energy 
( d ~ e p- We emphasize that e d — e p controls the corre- 
lation strength for charge-transfer materials in the ZSA 
framework, because the effective Coulomb repulsion in 
the one-band Hubbard model is controlled by e d — e p , 
not U dd [19]. 

To test our method, we use the extracted parame- 
ters for the canonical cuprate LSCO and explore the 
T = phase diagram as a function of doping. Our re- 
sults, shown in Fig. 4, are qualitatively similar to ex- 
periment. The calculations stabilize antiferromagnetism 
for low dopings x < 0.05, which gives way to a dome 
of d-wave superconductivity. The static order parameter 
A ~ (cc) T= o reaches a maximum A max near x ~ 0.12. 
We take the magnitude of A max as a proxy for the max- 
imum superconducting temperature T Cjmax . The zero- 
frequency limit of the anomalous self-energy E an is an 
additional indicator of superconductivity, which our re- 
sults show qualitatively follows the magnitude of the or- 
der parameter. 

We argue that although two independent low-energy 
parameters correlate with the experimental T c max , it is 
the charge-transfer energy that controls the variation in 
A max , and thus T Cimax , across the cuprate families. To 
address this issue, we take the most correlated cuprate, 
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FIG. 4. Calculated doping dependence for LSCO of the stag- 
gered magnetization S z = | (nf — n±) (diamonds) and static 
d-wave superconducting order parameter A ~ (cc) T =o- Opti- 
mal superconductiviting strength A max is obtained for dop- 
ing x opt ~ 0.13. The real part of the anomalous self-energy 
Re£ a "(aj = 0) (squares) follows qualitatively the order pa- 
rameter A. The calculations were performed at T — 30 K 
with c-DMFT and an ED impurity solver, using an 8-site dis- 
cretization of the bath. 



LSCO, and compute A max as we either (i) decrease e d —e p 
or (ii) increase t pp i . Fig. 5a shows that reducing the cor- 
relation strength for fixed t pp i enhances the order param- 
eter A, in agreement with the empirical trend in Fig. 3a. 
However, Fig. 5b shows that increasing t pp / across the 
physical parameter regime hardly modifies A max , in con- 
trast with the empirical trend in Fig. 3b. Further in- 
creasing t pp i to larger, unphysical values strongly sup- 
presses Tcmax- Thus, our calculations support the hy- 
pothesis that a larger hopping range r suppresses T Cjmax , 
in agreement with calculations on the one-band [16, 20] 
and three-band [14] models. 

The dependence of T c-max on the two controlled param- 
eters can be simply rationalized. For e d — e p , its large 
value in the strong correlation limit suppresses charge- 
fluctuations, rendering the residual superexchange inter- 
action between the doped holes weak, resulting in low 
superconducting temperatures. As we decrease e d — e p , 
superconducting tendencies increase as we pass through 
the intermediate correlation regime, until we reach the 
weak correlation limit. Although the ground state of 
the 3-band model for large U dd and e d — e p ~ has 
not been rigorously established, we expect the large ki- 
netic energy to suppress the effective interactions and 
thus superconductivity. Thus, we believe intermediate 
correlation strengths, a regime intimately related to the 
charge-transfer metal-to-insulator transition, is a crucial 
ingrediate for cuprate superconductivity. Turning to t pp i , 
we find that increasing this hopping amplitude lowers the 
van Hove singularity at (0, it) away from the Fermi level. 
The resulting decrease in density of states suppresses T c , 
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FIG. 5. Optimal superconducting order parameter A max of 
LSCO as we (a) decrease the charge-transfer energy ed — e v 
and (b) increase oxygen-oxygen hopping t PP i . Shaded are the 
physical ranges spanned by the cuprate families. 



an effect which simple methods capture [21]. 

In conclusion, we have used electronic structure meth- 
ods to identify the dependence of T c max on two funda- 
mental parameters: the charge-transfer energy ed — £p 
and inter-cell oxygen-oxygen hopping t pp >. We find that 
the position of the apical oxygen tunes both parameters, 
but the strength of superconductivity, A max is mainly 
sensitive to ed — e p . We expect future refinements to 
explain the remaining variability in T c max . Our work 
provides a natural interpretation of experiments where 
epitaxial compression in LSCO resulted in a remarkable 
enhancement of T c [22]. Epitaxy increases (icS-O ano - 
thus reduces ed — e p . Furthermore, our result provides 
microscopic insight into the multi-layer cuprates, such as 
Bi-2223: in addition to layer-dependent doping [23], the 
smaller value of the charge-transfer energy in the outer 
layers may explain the enhancement of superconductivity 
in the outer layers. It has been suggested theoretically 
and demonstrated experimentally [24] that proximity to 
a metallic layer reduces the charge-transfer energy. Us- 
ing this principle in heterostructure design should result 
in even higher transition temperatures. 
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SUPPLEMENTARY MATERIAL: TABLE OF PARAMETERS 

We summarize in Table I the parameters extracted via downfolding for the three-band model and discuss the details 
of the downfolding procedure. 

The charge-transfer energy td — e p is a localized, atomic-like quantity. Inherent in the downfolding procedure is a 
tradeoff between atomic character versus faithful representation low-energy bands. In order to preserve as much as 
possible the atomic character, we implemented the first step of the downfolding procedure described in Ref. 2. We 
chose as initial orbitals g n (r) the LDA+DMFT basis constructed in the manner described in Ref. 3. The downfolding 
procedure is robust: we cross-checked our results by using Wien2Wannier [4] and Wannicr90 [5] to perform the same 
downfolding procedure. Our code differs slightly in the choice of radial dependence of the trial orbitals g n {r). Again, 
in order to preserve the atomic character, we disabled the minimization of both spread functionals and did not use 
an inner window to constrain the Fermi surface. We find the extracted parameters differ by less than 5%. 

In order to connect with prior work [1], we compute the range parameter r ~ t'/t using Lowdin downfolding. 



Compound 


td - e p (eV) 


t pd (eV) 


tpp (eV) 


tpp' 


(eV) 


t'/t 


layers 


Co (A) 


Tc (K) 


(1) 


La2Cu04 


2.61 


1.39 


0.640 





103 


0.070 


1 


2.3932 


38 


(2) 


Pb 2 Sr 2 YCu 3 8 


2.32 


1.30 


0.673 





160 


0.108 


2 


2.3104 


70 


(3) 


C a2 Cu0 2 Cl2 


2.21 


1.27 


0.623 





132 


0.085 


1 


2.7539 


26 


(4) 


La 2 CaCu2C>6 


2.20 


1.31 


0.644 





152 


0.120 


2 


2.2402 


45 


(5) 


Sr 2 Nd 2 NbCu 2 Oi 


2.10 


1.25 


0.612 





144 


0.110 


2 


2.0450 


28 


(6) 


Bi 2 Sr 2 CuO e 


2.06 


1.36 


0.677 





153 


0.105 


1 


2.5885 


24 


(7) 


YBa 2 Cu 3 7 


2.05 


1.28 


0.673 





150 


0.110 


2 


2.0936 


93 


(8) 


HgBa 2 CaCu 2 6 


1.93 


1.28 


0.663 





187 


0.133 


2 


2.8053 


127 


(9) 


HgBa 2 Cu0 4 


1.93 


1.25 


0.649 





161 


0.122 


1 


2.7891 


90 


(10) 


Sr 2 Cu0 2 Cl 2 


1.87 


1.15 


0.590 





140 


0.108 


1 


2.8585 


30 


(11a) 


HgBa 2 Ca 2 Cu308 (outer) 


1.87 


1.29 


0.674 





184 


0.141 


3 


2.7477 


135 


(lib) 


HgBa 2 Ca 2 Cu308 (inner) 


1.94 


1.29 


0.656 





167 


0.124 


3 


2.7477 


135 


(12) 


Tl 2 Ba 2 Cu0 6 


1.79 


1.27 


0.630 





150 


0.121 


1 


2.7143 


90 


(13) 


LaBa 2 Cu3C>7 


1.77 


1.13 


0.620 





188 


0.144 


2 


2.2278 


79 


(14) 


Bi 2 Sr 2 CaCu 2 8 


1.64 


1.34 


0.647 





133 


0.106 


2 


2.0033 


95 


(15) 


Tl 2 Ba 2 CaCu 2 8 


1.27 


1.29 


0.638 





140 


0.131 


2 


2.0601 


110 


(16a) 


Bi 2 Sr 2 Ca 2 Cu30io (outer) 


1.24 


1.32 


0.617 





159 


0.138 


3 


1.7721 


108 


(16a) Bi2Sr2Ca2Cu30io (inner) 


2.24 


1.32 


0.678 





198 


0.121 


3 


1.7721 


108 



TABLE I: Tight-binding parameters of the three-band p-d model, containing the in-plane d x 2_ y 2 and P( x , y ) orbitals, for the 
hole-doped cuprates. The table is sorted by decreasing charge-transfer energies — e p . We have included the two nearest- 
neighbor (intra-cell) hoppings t p d and t pp as well as the inter-cell oxygen-oxygen hopping t pp i . Using the Lowdin procedure, we 
have integrated out the oxygen bands to arrive at a one-band model, from which we have extracted the ratio t'/t corresponding 
to the range parameter. We also include the distance between the in-plane copper and the apical atom eicu-o- For the bilayer 
and trilayer compounds, we display the distance to the apical oxygens from the Cu atoms in the outer planes. The last column 
displays the maximum transition temperature T™ ax for the corresponding optimally-doped compound. See Ref. 1 and citations 
therein for references to experimental work on structural determination and transition temperatures of the various cuprate 
families. 
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Beginning with the three-band model, 



/ 



H 



2t pd sin ~* 

y~2t pd sin -Jt —it pp sin ^ sin e p + 2t pp i cos fe y 



2t pd sin 



-2t pd sin -f 



I 2^-^^^ cos 



4t pp sin ^ sin ^ 



we integrate out the oxygen bands to arrive at the effective one-band hamiltonian, 



H cS (uj) -e d + t pd ■ —2 , 

2^i=o Bi(u))ai(\t) 



where we have defined the Fourier harmonics as 



a (k) 
ai(k) 
a 2 (k) 



= 1 



— 2(C0S fcr + COS fey) 

4 cos fer cos k„ 



and the coefficients are 

-4o = 4t pd (uj - e p + 2t pp ) B = (u- e p ) 2 - t 2 pp 

Ai = t pd (uj — e p + 2t pp i + 4t pp ) B\ = t pp i (ui — e p ) — 2t pp 

A-2 = 2t pd (t pp + t pp ') B 2 = tpp, - t p p. 

Taking advantage of the fact that B\ / Bq and B2 /Bo are small, we expand out the denominator and collect coefficients 
of the Fourier harmonics to arrive at the range parameter 

I = B A 2 - A B 2 
r ~ * BoA, - A B 1 w=£f 1 > 

This procedure preserves the Fermi surface and faithfully represents the low-energy band-structure. Since we priori- 
tized the faithful representation of the atomic quantities over the non-local hopping parameters, our values of r ~ t' /t 
are smaller than those found in Ref. [1]. However, the trends remain unchanged. 



SUPPLEMENTARY MATERIAL: NUMERICAL METHOD 



We solve the following three-band hamiltonian using c-DMFT: 

H = ^2 t"f<$ a<T Cjp„ + ^2 e a n ioia + U dd ^ n ld ^n ld ^. (2) 

iaj/3(7 iaa ia 

We use in this work the realistic set of parameters shown in Table (I) . We choose a 2 x 2 cluster in the cellular form [6] . 
Cluster DMFT improves on the single site DMFT by adding a non-local self-energy. In Fig. l.a we show the 2x2 
copper plaquette used as a unit cell through the calculations. The lattice Green's function of the four Copper site 
plaqucttc is given by: 

Gk(iWn) = (iu) n + A* - H k - E(iw n )) _ , (3) 

where Hk is the Fourier transform of the uncorrelated part of the Hamiltonian defined in Eq. (2). S is the cluster 
self-energy matrix being nonzero only for the matrix elements connecting the d x 2_ y 2 orbitals. 

The self energy matrix in Eq. (3) is obtained by solving an 2 x 2 impurity Anderson model (shown in Fig. l.b) 
subject to the DMFT self-consistency condition: 

(iu> - E lmp - S(zc^) - AM) = P I -±r G kM ] > ( 4 ) 
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FIG. 1: (Colors online) a) The d-p theoretical lattice model contains the p< x ,y) orbitals of the in-plane oxygen atoms (small 
circles) and the d x 2_ y 2 orbital of the copper atoms (large circles). The figure shows the four copper plaquette unit-cell used 
in the cellular DMFT calculations, b) The associated Anderson impurity model (AIM) contains four impurities (circles), each 
of them is independently connected to a bath discretized in eight sites (stars). There is not direct hybridization between the 
impurities, but the latter are connected by second order process through the bath (dashed lines). The sites of the bath are 
connected by direct and long-range hoppings. 



where the sum runs over the reduced Brillouin Zone (BZ), and P is projecting the averaged green function onto the 
impurity cluster subspace. 

In this work we use the exact diagonalization impurity solver algorithm [7]. To solve the cluster impurity problem, 
we express it in the form of a Hamiltonian Hi mp with a discrete number of bath orbitals coupled to the cluster and 
use the Lanczos algorithm to converge the ground state of the Hamiltonian and the lower states of the spectrum. The 
ED method in conjunction with c-DMFT has been widely used for the one-band model [8-10]. 

The Anderson Impurity Model (AIM) is defined by : 



The fermionic operators a\ nn (a mn ) creates (destroys) a particle in the bath, and the fermionic operators c\ nn (c mn ) 
creates (destroys) a particle in the cluster of impurities. The indices m, n are running over the bath sites, and the 
index i is running over the impurity sites, the sites of the bath are connected by long-range hopping matrix elements 
through the particle-hole (particle-particle) channel e™ (e a ), the non-correlated sites of the bath are also connected 
to the correlated impurities by the matrix elements V m i, the onsitc repulsion at the impurity sites is U (equal to the 
Coulomb repulsion of the Copper site Ud), and \i is the chemical potential. We define e as the extended matrix which 
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contains the normal e n and anomalous e a blocks in the Nambu basis: 



e = 



( 





a 



) 



(0) 



The Weiss field Q(iu> n ) = ito n — A(iu) n ) — E imp is constructed from the parameters of the AIM: 



A(tw n ) = V f (i 




(7) 



The parameters of (5) are determined by imposing the self-consistency condition in Eq. (4) using a conjugate gradient 
minimization algorithm: 



Where a/3 are running over the matrix elements, ojq = 20 is a hard cutoff on the summation and A ED is the function 
(7) of the Hamiltonian parameters. The fitting procedure is not exact due to the discretization of the bath and is an 
additional approximation to the c-DMFT scheme. In this work we considered a bath discretised with 8 energy levels. 
Finally, once the Hamiltonian parameters are obtained by the fitting procedure, we obtain the low energy spectrum 
by the Lanczos procedure. We impose an energy cutoff E max such that the Boltzman weig ht e -ft E ™*- E d < 0.001, 
where Eq is the ground state energy. We discard all the eigenstates which have an energy larger than the cutoff 
Ei — E > E rnax . Once the eigenstates are obtained we compute the Boltzman weighted average to get the dynamical 
and static observables. 

In this work we consider two different instabilities: i) the superconducting phase, and ii) the long-range magnetic 
ordered phase. The former is computed in the Nambu basis, and the Hilbert space is block diagonalized by the spin 
S z quantum numbers, and the latter is obtained in the tensor product of the up and down spins. Since the number 
of particles is not a good quantum number in the superconducting phase, we work at fixed chemical potential. For 
the magnetic phase we found a better convergence when working at fixed density, with a free chemical potential. 

For the determination of the phase diagram we used physical observables readily available from the 2x2 cluster 
of impurities, such as the staggered magnetization S z — \ (n-j- — 714.), the superconducting order parameter A = 
(( c 1 c 2))(t=o) (where 1 and 2 are nearest neighbor links of the impurity plaquette) and the anomalous self-energy at 
zero frequency S an = = 0). We emphasize that the computed order parameters do not rely on any additional 

procedure, such as the E-periodization [6], which interpolates and extrapolates the discrete cluster quantities to the 
continuum in k-space. 
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